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I .  INTRODUCTION 


An  antiarraor  munition  is  being  developed  and  tested  by  the  US  Army  Arma¬ 
ment  Research  and  Development  Command  (ARRADCOM) .  The  lead  laboratory  for  the 
program  is  ARRADCOM' s  Large  Caliber  Weapon  System  Laboratory  (LCWSL) .  The 
system's  analysis  and  sensor  technology  research  is  ongoing  at  the  Ballistic 
Research  Laboratory  (BRL) ,  Aberdeen  Proving  Ground,  MD.  As  part  of  this  BRL 

effort  the  millimeter  wave  (MMW)  radiometric  sensor  was  tested  in  Aigust  1978 

at  the  North  East  Test  Site  of  the  Rome  Air  Development  Center  in  Rome,  NY. 

This  electromagnetic  radiation  data  collected  in  Rome,  NY,  is  related  to 
the  absolute  temperature  of  the  ground  surface  at  the  sensor's  focus  point,  and 
is  referred  to  as  radiance  n-  The  functional  relationship  between  absolute 
temperature  and  the  radiance  is  given  in  equation  (1) : 

n  =  2kT/x2  Cl) 

where  k  is  Boltzmann's  constant 

T  is  the  absolute  temperature 
and  X  is  the  wavelength. 

The  radiometric  data  were  collected  over  field  and  wooded  areas,  at  five 
different  slant  ranges  from  30  to  150  meters  (see  Figure  1) .  The  MMW  sensor 
was  mounted  on  a  helicopter  at  a  squint  angle  of  30?,  measured  from  perpen¬ 
dicular  to  the  ground.  The  helicopter  was  then  flown  at  a  ground  speed  of 
60  knots.  The  MMW  sensor  had  a  spin  rate  of  4  revolutions  per  second  and 
recorded  data  continuously.  Then,  the  data  was  digitized  by  sampling  approxi¬ 
mately  two  thousand  equally  spaced  observations  per  second. 

Representative  plots  of  the  August  1978  MMW  data  are  presented  in 
Figure  2  and  Figure  3.  Figure  2  is  a  typical  terrain  (non- target)  response, 
whereas  Figure  3  is  a  response  with  a  target.  The  absolute  temperature  T  in 
Equation  (1)  has  been  adjusted  in  Figure  2  and  Figure  3  by  subtracting  the 
average  ground  temperature  plus  an  arbitrary  offset.  Hereafter,  the  symbol 
T  will  be  referred  to  as  the  contrast  temperature.  The  analysis  of  the 
system's  response  to  a  target  (Figure  3)  will  be  discussed  in  Section  IV. 

The  purpose  of  the  sensor  is  to  detect  the  target  and  then  aim  the 
weapon  system  under  battlefield  conditions.  Hence,  the  response  of  the 
MMW  sensor  to  varying  background  conditions  will  affect  the  weapon  system's 
ability  to  detect  targets. 

The  radiometric  non-target  and  target  observations  were  modeled  using 
the  Box  and  Jenkins  approach.  The  analysis  demonstrates  the  ability  of  the 
ARIMA  model  to  characterize  the  radiometric  data. 
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Figure  1.  NORTH  EAST  TEST  SITE  OF  THE  ROHE  AIR  DEVELOPMENT  CENTER  IN  ROME,  N.Y. 
AUGUST  1978 
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CONTRAST  TEMP  vs.  TIME 


II.  THE  BOX  AND  JENKINS  TIME  SERIES  APPROACH 

The  Box  and  Jenkins1  Autoregressive  Integrated  Moving  Average 
(ARIMA  (p,d,q))  model  is  used  to  characterize  many  types  of- business, 
economics  and  engineering  observations.  A  requirement  to, develop  a  model  of 
the  MMW  sensors'  responses  to  backgrounds  (terrain)  has  existed  at  BRL 
for  some  time.  A  new  approach  using  the  Box  and  Jenkins  time  series  metho¬ 
dology  was  used  to  model  the  sensors'  responses. 

The  ARIMA  (p,d,q)  model  is  presented  below: 

cj)p(B)(l-B)d(Zt-y)  =  0  (B)at  (2) 

where  B  is  the  backshift  operator  such  that  BZ^  =  Z 

^p(B)  is  a  polynomial  in  B  of  order  p  and  <Jk  are  autoregressive  parameters, 

Qq(B)  is  a  polynomial  in  B  of  order  q  and  9^  are  moving  average  parameters, 

p,d,q  are  non-negative  integers,  and  refer  to  the  order  of  the  respective 
autoregressive  operators,  the  difference  and  the  moving  average 
operator, 

at  are  random  shocks  (white  noise)  assumed  to  be  independently 

2 

distributed  normal  variates,  N(0,cr  )  at  time  t,  and  in  this  application 

Vi 

\  are  values  of  contrast  temperature  at  the  t  time  step. 

A  series  of  five  hundred  observations  were  analyzed  and  used  to 
estimate  the  Autocorrelation  Function  (ACF)  and  Partial  Autocorrelation 
Function  (PACF)  out  to  40  lags  or  time  steps  (see  Table  1  and  Figure  4) .  The 

estimated  mean  (y)  and  standard  deviation  (a)  of  this  500  sample  time  series 
are  -1 . 603  and  5, 712,  respectively. 


In  building  a  dynamic  time  series  model,  three  steps  are  required; 

(1)  tentatively  identify  the  model,  (2)  fit  the  data  to  the  model,  and’ 

(3)  perform  a  diagnostic  check  for  lack  of  fit.  The  identification  step 
requires  an  overall  view  of  the  data  structure.  In  this  case  a  damping 
sinusoidal  structure  for  the  ACF  and  two  significant  spikes  at  the  first 
and  second  PACF  (see  Figure  4)  implies  that  an  autoregressive  model  of  order 
p=2  be  initially  entertained.  Hence,  the  ARIMA  (2,0,0)  model  was  fitted. 


ARIMA  (2,0,0)  Initially  Entertained 


(1  -  ^B  -  <j)2B2)  (zt  -  y)  =  at 
where  the  estimated  parameters  are  y  =  -1.636 


4^  =  1.252 

✓s 

<f>2  =  -0.450. 


(3) 


Box,  E.  E.  P.  and  Jenkins,  G.  M.,  Time  Series  Analysis:  Forecasting  and 
Control,  San  Francisco,  CA,  Holden-Day,  1970.  - - 
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TABLE  1.  ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION  FUNCTIONS 

MMW  RADIOMETRIC  READING!  2000  observations  per  second 

Autocorrelations  (ACF) 


500  Observations 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

z 

Lags  1-10 

.86 

.63 

.46 

.35 

.27 

.21 

.17 

.15 

.15 

.15 

11-20 

.17 

.22 

.27 

.33 

.37 

.40 

.39 

.36 

.32 

.28 

21-30 

.23 

.17 

.11 

.07 

.04 

.04 

.06 

.12 

.18 

.23 

31-40 

.26 

.31 

.36 

.37 

.31 

.23 

.16 

.10 

.04 

-.03 

to 

Lags  1-10 

.34 

-.21 

-.22 

-.11 

-.08 

-.07 

-.08 

-.06 

-.02 

-.06 

11-20 

-.09 

-.04 

.02 

.01 

.09 

.12 

.06 

.03 

.04 

.01 

21-30 

.02 

.01 

-.03 

-.10 

-.06 

-.11 

-.11 

-.02 

.06 

.04 

31-40 

-.05 

.01 

.15 

.22  • 

.12 

-.07 

-.05 

.04 

-.00 

-.12 

Partial  Autocorrelations  (PACF) 

500  Observations 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Lags  1-10 

.86 

-.43 

.25 

-.09 

.03 

.01 

.02 

.04 

.03 

.02 

11-20 

.12 

.09 

.08 

.09 

.11 

-.02 

.03 

.03 

-.01 

-.02 

21-30 

-.00 

-.06 

-.03 

-.03 

.03 

-.05 

.11 

.04 

.02 

.01 

31-40 

.06 

.14 

.01 

-.06 

-.09 

-.04 

.06 

-.10 

-.06 

-.04 

Lags  1-10 

.34 

-.37 

.01 

-.12 

oo 

o 

• 

* 

-.09 

-.11 

-.08 

-.07 

-.15 

11-20 

-.12 

-.10 

-.11 

-.12 

.02 

-.05 

-.04 

-.00 

.01 

-.01 

21-30 

.05 

.01 

.02 

-.05 

.04 

-.14 

-.05 

-.04 

-.01 

-.08 

31-40 

-.14 

-.01 

.05 

.08 

.04 

-.08 

.10 

.04 

.03 

-.05 
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yz  -  -1.603 

A 

Oz  =  3.712 


I  .0  T 


0.5  « 


*1 


0.0 


-0.5  «■ 


-1 .0 


PACF(Zt) 


Li  1  1  1  J  t  .  .  f  1 


10 


1  j  «  *  l 

20 


30 


tj  *  j  n 

H0 


Figure  4.  Estimated  Autocorrelation  (rfc)  and  Partial  Autocorrelation 
Functions  (<J>k)  for  the  kth  Lag. 
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The  autocorrelation  function  of  the  residual,  a  =  Z  -  Z^.,  of  the 

ARIMA  (2,0,0)  model  was  checked  for  lack  of  fit.  These  residual  ACF  are 
listed  in  Table  2,  where  the  estimated  residual  mean  and  standard  deviation 

/s  A 

are  y  =  0.00175  and  a  =  1.681.  The  ACF  of  the  residuals  at  lags 
at  at 

k=l  and  k=2  indicated  some  remaining  residual  structure  (see  Figure  5) . 

The  cutoff  of  the  residual  ACF  suggests  a  possible  need  for  a  moving 
average  term  in  the  model . 

Hence,  the  ARIMA  (2,0,1)  model  was  entertained  in  order  to  remove  the 
spikes  in  the  residual  ACF.  The  ARIMA  (2,0,1)  model  is  as  follows: 


ARIMA  (2, 

0,1) 

(1  - 

-  <1>2B2)  (Zt  ■ 

-  y) 

=  (i  - 

9lB)at 

where  the  estimated 

parameters 

are 

y  = 

-1.622 

/\ 

♦i  = 

0.886 

*2  = 

-0.138 

II 

t— l 

<  © 

-0.506  . 

The  addition  of  the  moving  average  term,  0^,  removed  the  autocorrelation 
spikes  for  lags  k  =  1  and  k  =  2  (see  Table  3) .  Inspection  of  the  residuals 

A  /\ 

with  mean  u  =  0.0009  and  standard  deviation  a  =  1.616,  where  the  sample 
t  at 

ACF  and  PACF  values  are  small,  indicates  a  good  fit  (see  Figure  6) .  The  95% 

A 

confidence  interval  for  the  second  estimated  autoregressive  parameter,  <p^  =  -0.138, 

overlapped  the  zero  point  (-0.285,  0.0086).  This  suggested  the  possible  removal 
of  this  term  from  the  model.  Based  on  this  information  and  the  principle  of 
parsimony  in  the  use  of  parameters,  the  second  autoregressive  parameter,  <f>~, 
was  removed  and  the  ARIMA  (1,0,1)  model  was  then  considered.  2 

The  ARIMA  (1,0,1)  model  is  as  follows: 

ARIMA  (1,0,1) 

(1  -  <J>1B)(zt  -  y)  =  (1  -  61B)at  (5) 

-A. 

where  the  new  estimated  parameters  are  p  =  -1.618 

(f)1  =  0.7553 

A 

0X  =  -0.5960  . 
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TABLE  2 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION  FUNCTIONS  OF  RESIDUALS 

(ARIMA  (2,0,0)) 


Residual  Mean  ■  .001747 

Residual  Standard  Deviation  =  1.6809 


*k 


ACF 


500  Observations 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Lags  1-10 

.12 

-.21 

.01 

.10 

.07 

.04 

.01 

.01 

.08 

.02 

11-20 

-.02 

.04 

.09 

.01 

.11 

.14 

.05 

.04 

.08 

.04 

21-30 

.06 

.05 

.04 

-.05 

.06 

-.01 

-.05 

.04 

.10 

.08 

500  Observations 

1 

2 

3 

PACF 

4 

5 

6 

7 

8 

9 

10 

Lags  1-10 

.12 

-.22 

.07 

.05 

.07 

.06 

.02 

.02 

.07 

O 

O 

* 

11-20 

-.00 

.04 

.06 

t 

• 

o 

o 

.14 

.10 

.06 

.06 

.07 

.01 

21-30 

.05 

.01 

.03 

00 

o 

» 

1 

.05 

-.09 

-.04 

-.02 

.06 

.05 
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ACF  (a) 
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Figure  5.  Plot  of  the  Estimated  ACF  (rk)  and  PACF  (<J>k)  of  Residuals 
(ARIMA  (2,0,0))  for  the  kth  Lag. 
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TABLE  3 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION  FUNCTIONS  OF  RESIDUALS 


(ARIMA  (2,0,1)) 


Residual  Mean  =  0.0009 

Residual  Standard  Deviation  «=  1.6169 


500  Observations 

 12  3 


Lags  1-10 

.00 

-.02 

-.03 

11-20 

-.02 

.01 

.09 

21-30 

.07 

.01 

.05 

500  Observations 

1 

2 

3 

Lags  1-10 

.00 

-.02 

-.03 

11-20 

<N 

O 

• 

1 

.02 

• 

o 

00 

21-30 

.07 

.01 

.06 

ACF 


i  5_ 6  7  8  9  10 


05 

.01 

.02 

-.01 

-.01 

.05 

o 

o 

• 

1 

00 

.11 

.11 

.06 

.04 

.09 

.02 

09 

.05 

-.05 

-.04 

.02 

.06 

.07 

PACF 

4  5 

6 

7 

8 

9 

10 

05 

.01 

.02 

o 

p 

• 

1 

-.02 

.05 

-.01 

00 

.12 

.12 

.06 

.05 

.09 

.02 

09 

.04 

-.06 

-.06 

-.01 

.03 

.05 

19 


1  -0  T 


.5  «■ 


0.0 


—0.5  «■ 


-1 .0  A 


PACF  (a) 


-i— «.  ji  .  .  y 


■  *  ii  * .  i .  i 


J-T-*  ...»  1 


IB 


20 


"T~TT 


30 


Figure  6.  Plot  of  the  Estimated  ACF  (rk)  and  PACF  (<f>k)  of  Residuals 
(ARIMA  (2,0,1))  for  the  kth  Lag. 
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Both  the  ACF  and  PACF  of  the  ARIMA  (1,0,1)  residuals  indicated  an 
excellent  fit.  Table  4  and  Figure  7  show  that  the  residual, 

fat  =  zt  -  zt>,  contains  no  further  structure.  Additional  attempts  to 

better  characterize  this  time  series  were  unsuccessful. 

A  summary  of  the  models  and  their  estimated  parameters  is  presented 
in  Table  5. 


III.  APPLICATION  OF  THE  ARIMA  (1,0,1)  MODEL  TO  ADDITIONAL  NON-TARGET  MMW  DATA 

In  Section  2,  the  rationale  for  the  selection  of  the  ARIMA  (1,0,1)  model 
was  presented.  In  this  section,  additional  non-target  cases  are  analyzed. 

The  ARIMA  (1,0,1)  model  was  found  to  be  adequate  for  modeling  these  additional 
sets  of  non-target  MMW  data  collected  in  Rome,  N.Y. 

In  Table  6,  both  the  ACF  and  PACF  of  additional  samples  of  500  observa¬ 
tions  at  different  distances  (30,  60,  90,  120  and  150  meters)  were  estimated. 
Plots  of  these  ACF's  and  PACF's  are  presented  in  Figure  8.  The  similarity  in 
the  correlograms  raised  the  possibility  of  modeling  all  the  non-target  data 
with  the  same  ARIMA  model.  Hence,  Table  7  was  generated  from  30  additional 
sets  of  data  using  the  ARIMA  (1,0,1)  model. 

/\ 

This  analysis  suggested  that  the  estimated  means  (y)  are  not  fixed,  ^ 
whereas  the  autoregressive  parameter  (<f>.)  and  the  moving  average  parameter  (0.) 

varied  only  slightly.  The  white  noise  estimated  parameters  (y  ,  cr  )  were  consis- 

A,  3.  cl 

tently  similar  for  those  cases  investigated.  That  is,  y  0.0  and 

A  cl 

°a  **  C1-62  t0  i*79)*  further  suggesting  the  utility  of  this  noise  model 
approach.  A  closer  look  at  the  residuals,  at,  indicated  a  lack  of  any  consis¬ 
tent  pattern  after  fitting  the  ARIMA  (1,0,1)  model.  Figure  9  shows  this  lack 
of  structure  in  the  residuals,  indicating  again  the  ability  of  the  ARIMA  (1,0,1) 
model  to  characterize  the  balance  of  the  1978  MMW  non-target  data.  A  white 
noise  model  was  used  to  simulate  the  sensor's  characteristics  for  the  different 
cases.  The  ARIMA  model  used  for  this  purpose  is  that  of  Equation  (6). 

A  A  A  A 

zt  =  (1  -  4>)y  +  at  -  6at_1  +  <pzt_1  (6) 


^2 

where  a  =  N(0,  o  )  and 

L  cl 

AAA 

(y,<J>,0)  are  the  estimated  parameters. 

A  plot  of  one  such  simulated  case  is  shown  in  Figure  10  for  a  comparison  with 
the  actual  data  plotted  in  Figure  2. 
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TABLE  4 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION  FUNCTIONS  OF  RESIDUALS 

(  ARIMA  (1,0,1)) 


Residual  Mean  *  0.00039 

Residual  Standard  Deviation  *  1.623 
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Figure  7.  Estimated  ACF  (r^)  and  PACF  ((j)^)  of  Residuals 
(ARIMA  (1,0,1))  for  the  kth  Lag. 
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TABLE  6 


AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION  FUNCTIONS  OF  Z. 

«  t 


500  Observations 


ACF 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Zt  Lags  1-10 

.86 

.63 

.46 

.35 

.27 

.21 

.17 

.15 

.15 

.15 

(150  m)  11-20 
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.40 
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.39 

.36 

.32 

Zt  Lags  1-10 

.86 

.64 
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.40 

.35. 

.30 

.24 

.19 

.18 

.21 

(120  m)  11-20 

.27 
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.37 
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.50 
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.42 

.36 

Zt  Lags  1-10 

.85 

.60 

.44 

.35 

.28 

.20 

.15 

.13 

.14 

.18 

(90  m)  11-20 

.21 

.25 

.32 

.41 

.48 

.53 

.52 

.46 

.38 

.30 

Z^  Lags  1-10 

.87 
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.31 
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.07 
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Zt  Lags  1-10 

.86 
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Figure  8.  The  ACF  and  PACF  of  (Zt) 


TABLE  7 

A  SUMMARY  OF  ADDITIONAL  MMW  DATA  FITTED  TO  THE  ARIMA  (1,0,1)  MODEL 


August  2,  1978  (Field) 

A 

♦l 

.7329 
.96376* 

.76154 
.79747 

.72707 
.68369 

.75301 
.7572 

.7552 
.68908 
.66975 

.732691 
.040613 

.753697) 
.079614) 


August  4,  1978  (Field) 


30  Meters 

(a) 

-  9.3985 

.89142 

(b) 

-  7.8494 

.79993 

60  Meters 

(a) 

6.8121 

.87228 

(b) 

8.4699 

.  805B9 

90  Meters 

(a) 

23.055 

.88165 

(b) 

23.037 

.86254 

120  Meters 

(*) 

1.2813 

.91947 

(b) 

1.4580 

.92505 

150  Meters 

(a) 

-14.029 

.83538 

(b) 

-14.726 

.87666 

Mean 

M/A 

.867027 

St 

.  Dev. 

N/A 

.042567 

August  4, 

1978  (Trees) 

30  Meters 

(a) 

-  1.2039 

.73473 

(b) 

-  1.4959 

.81564 

60  Meters 

(a) 

-  2.6053 

.82087 

(b) 

-  2.0881 

.75270 

90  Meters 

(a) 

16.155 

.72748 

(b) 

15.704 

.83907 

120  Meters 

(a) 

-  7.6494 

.78241 

(b) 

-  7.2249 

.81785 

150  Meters 

(a) 

8.5403 

.85605 

(b) 

9.0061 

.79945 

Mean 

N/A 

.794625 

St 

.  Dev. 

N/A 

.044026 

ARJMA  (1.0,1) 

(W,B)(Zt 

-Vz)  •  (1- 

e,B).t 

- 

A  • 

A 

fll 

°a 

*•(27 ) 

-.67395 

-.00078 

1.6542 

49.21 

*  -.55912 

-.000895 

1.6323 

56.06 

-.54864 

-.000674 

1.5790 

58.89 

-.60273 

-.000316 

1.6816 

47.83 

-.63197 

.000253 

1.6216 

65.50 

-.67054 

-.0003039 

1.6275 

56.78 

-.56238 

.0001229 

1.5556 

70.70 

-.57466 

-.0000809 

1.5959 

54.47 

-.5959 

.000399 

1.6231 

45.80 

-.61738 

.000671 

1.6515 

42.67 

-*  67832 

.0001374 

1.5796 

74.36 

-.610632 

-.0001333 

1.618355 

t  .047867 

.000507 

0.037667 

- 

i ) 

ARIHA  (1,0.1) 

-.57673 

-.0002821 

1.7345 

42.50 

-.61042 

-.0000485 

1.7487 

32.29 

-.56256 

-.0005989 

1.7910 

45.25 

-.66890 

.0002418 

1  .8386 

47.28 

-.82222 

-.0000369 

1  .8039 

55.20 

-.56396 

.00031859 

1.7622 

51.01 

-.53543 

.0003456 

1.7566 

55.01 

-.62257 

.0000953 

1.7982 

38.34 

-.67237 

-.000327 

1.7767 

39.40 

-.62128 

-.0004968 

1.9116 

55.95 

-.605644 

- . 00007B89 

1.7922 

.045512 

.00033697 

0.051868 

- 

ARIMA  {1,0, 

1) 

-.63485 

-.0001325 

1.6602 

69.56 

-.59693 

-.0000552 

1.747 

47.10 

-.59368 

-.0004588 

1.7232 

37.90 

-.63814 

-.0001127 

1.7423 

58.36 

-.68935 

.0002787 

1.7499 

42.76 

-.63629 

-.0004865 

1.8845 

74.76 

-.62126 

.0001219 

1.8704 

45.77 

-.62218 

.000946 

1.8263 

70.69 

-.59373 

-.0002017 

1.8195 

53.10 

-.61361 

-.001264 

1.8421 

77.85 

-.624002 
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1.78654 

.028726 

.00057 

0.07243 

30  Meters 

(a) 

-20.43 

(b) 

-24.12 
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(a) 

-17.185 

(b) 

-17.456 

90  Meters 

(a) 

-  7.5447 

(b) 

-  9.2541 
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(a) 

-  7.514 

(b) 

-  7.909 

150  Meters 

(a) 

-  1.618 

(b) 

-  1.7179 

(c) 

-  0.B284 

He.n 

N/A 

St.  Dev. 

N/A 

27 


1.0 


ACFC150  M) 
2  Aug  78 

-1.0 

1.0 


ACF(120  M) 
2  Aug  78 


.....it,  ilt 

I  1  * 


-1.0 

1.0 


ACF(90  M)  .  ,  .  .  .  .  .  I  1  i 

2  Aug  78  1 

-1.0  ■■ 


1.0 


ACF(60  M) 
2  Aug  78 


L-* - 1  ,  t  .  *  i-l—M—l  -i—i 


TT 


x 


-1.0 

1.0 


ACF  (30  M)  i  .  ,  -  i _ l  -l-l.  .  .  « 

2  Aug  78 

-1.0  .. 


Figure  9.  The  Autocorrelation  Function  of  ARIMA  (1,0,1)  Residuals 
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Figure  10.  Modeled  MMW  Time  Series 


In  summary,  the  MMW  non-target  data  collected  in  August  1978  was 
analyzed  using  the  Box  and  Jenkins  Time  Series  approach.  This  analysis  indi¬ 
cated  an  ARIMA  (1,0,1)  model.  This  particular  ARIMA  model  characterizes  the 
data  remarkably  well.  What  is  more  remarkable  is  the  consistent  behavior 

of  both  the  estimated  parameters  C<J>1,  0j)  and  the  white  noise  parameters 

(ya»  °a)  •  is  recommended  that  the  ARIMA  (1,0,1)  model  be  used  to  describe 
the  non-target  data. 


IV.  THE  USE  OF  INTERVENTION  ANALYSIS  TO  CHARACTERIZE  TARGETS2 

In  the  previous  section,  an  ARIMA  model  was  used  to  characterize  the  non 
target  data.  An  additional  requirement  is  to  model  the  data  when  a  target 
is  present.  The  approach  used  is  referred  to  as  intervention  analysis.  By 
intervention  we  mean  that  a  target  intervenes  or  projects  itself  on  the  non¬ 
target  data  causing  changes  to  the  time  series. 

In  Figure  3  a  typical  plot  of  the  MMW  sensor's  responses  to  a  target 
is  shown.  This  target  response  is  denoted  by  the  negative  dip  in  the 
stationary  time  series  at  approximately  t  =  275.  In  this  section,  a  dynamic 
intervention  model  is  being  proposed  to  describe  the  targets'  responses 
This  dynamic  model  is  illustrated  in  Figure  11.  For  the  present,  the  target 
will  be  denoted  by  a  pulse  indicator  function. 


P 


t 


1, 

0, 


Targets 

Otherwise 


(7) 


where  T  is  the  number  of  pulses  in  the  signal.  That  is,  when  the  sensor  detects 
the  occurrence  or  nonoccurrence  of  a  target  the  variable  takes  on  the  value 
1  or  0.  The  input  variable  is  denoted  by  Ct  =  WoPt(T).  This  dynamic  model 

takes  the  form,  Zt  =  +  N't,  where  the  Yt  represents  the  additional 

effect  of  an  intervention  over  the  stochastic  noise  (Nt) . 

In  the  attempt  to  model  targets,  an  intervention  dynamic  model  of  the 
following  form  is  proposed: 


U  -  WjB)  (l  _  0  R) 

Zt  =  (1  -  SjB)  Wo  PtCT)  +  IT  -  ^1)  at 


Where  P  (T) 


■{ 


1,  Targets 
0,  Otherwise 


and 


>t  =  N(0-  O 


C8) 


Box,  G.  E.  P. 
Economic  and 
Association, 


and  Tiao,  G.  C.,  "Intervention  Analysis  with  Applications  to 
Environmental  Problems,"  Journal  of  the  American  Statistical 
Vol.  70,  No.  349,  March  1975. 
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Figure  11.  Intervention  Model  for  a  Dynamic  System  with  Superimposed  ARIMA  Model 


This  structure  is  being  proposed  because: 

(1)  In  Sections  II  and  III  the  noise  is  modeled  by  (1  -  <j>  B)N  =  (1  -  0,83a 

and  the  estimated  parameters  are  approximately  it  It 

A 

<{)  %  .80 

/S 

0  »  -  .62 

/N 

Pa  «  0.00 

/\ 

a  «  1.78 
a 

(2)  We  take  the  approach  that  the  raw  signals  of  both  background  noise 
and  target  will  filter  through  the  identical  MMW  sensor.  Hence,  both  signals 
are  likely  to  have  similar  dynamic  structure.  The  above  information 
resulted  in  the  following  time  series  intervention  model  where  the  previous 
estimated  parameters  are  substituted  into  Equation  (8) . 


(1  +  0.62B)  „ 
Cl  -  0.8B)  o 


PtCT) 


+  y  + 


(1  +  0.62B) 

Cl  -  0.8B3  at 


where 


Target 

otherwise 


W  <  0. 
o 


(9) 


Before  continuing,  a  look  at  the  dynamic  response  to  a  pulse  input 
under  the  proposed  intervention  structure  would  be  helpful.  (See  Table  8.) 


Let 


(1  -  WXB) 

Yt  =  '  (1  -  sJb) 


W  (T) 
o  tv  * 


(10) 


or  Yt  =  6iYt_i  +  W0Ptm  "  WlWoPt-l(T)  CH) 

The  initial  responses  will  be  a  step  of  size  W  ,  followed  by  (T-l) 

additional  increasing  steps.  This  dynamic  intervention  model  would  be  at  its 
maximum  value  at  the  Tth  step  (the  pulse  length)  and  would  start  to  decay  at 
a  rate  <5^  If  Wq  were  negative  (W  <  0)  the  response  of  the  pulse  input 

would  also  be  negative  similar  to  the  actual  targets'  responses.  That  is,  the 
response  would  occur  in  the  negative  direction  instead  of  the  positive  direction 
as  shown  in  Table  8.  The  number  of  pulses  needed  to  generate  these  dynamic 
responses  can  almost  be  interpreted  as  the  number  of  times  the  MMW  sensor 
detects  the  presence  of  the  intervention  (target).  The  size  of  the  adjacent 
response  can  be  thought  of  as  a  function  of  W  ,  the  intensity  of  the  given 
target.  That  is,  if  the  adjacent  step  responses  are  large,  the  target  tempera¬ 
ture  characteristics  are  high.  Hence,  there  exists  good  rationale  for  the 
parameters  Wq  and  the  dynamic  intervention  model,  which  are  summarized  below: 
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TABLE  8.  THE  DYNAMIC  RESPONSES  OF  A  PULSE  INPUT  FOR  THE 
PROPOSED  INTERVENTION  STRUCTURE. 


Yt  "  Vt-l  +  W0  Pt<T)  -  WTW0 


where  6-j  =  0.8 
WQ  =  1.0 

W1  =  -0.62 


INPUT 


RESPONSES 


T  =  1 


Uni. 


T  =  3 


III 
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TABLE  8. 


(Continued)  THE  DYNAMIC  RESPONSES  OF  A  PULSE  INPUT 
FOR  THE  PROPOSED  INTERVENTION  STRUCTURE. 


INPUT 


RESPONSES 


T  =  5 


T  =  6 


(i)  The  presence  of  a  target  is  indicated  by  a  pulse  (intervention) 
of  length  equal  to  that  of  a  target  signal. 

(ii)  The  raw  target  signal,  like  the  background  noise, 

should  behave  in  a  similar  dynamic  structure  (same  ARIMA  model) 
(1  -  61B)Zt  =  (1  -  WjBJCJ.  ' 

(iii)  The  number  of  pulse  inputs  is  an  indication  of  how  long  the 
target  is  being  sensed  by  the  MMW  sensor. 

(iv)  The  step  "Wo"  is  an  indication  of  the  contrast  temperature 
T(°K)  which  is  a  characteristic  of  the  target  (intervention). 

This  representation  is  in  its  formative  stage  and  improvements  are  planned  as 
more  target  data  become  available.  It  is  believed  that  =  WQP  (T)  is  not 

in  reality  a  constant,  but  a  function  of  the  size  and  temperature  characteristics 
of  the  intervention  (target) . 

This  intervention  model  was  used  to  simulate  a  number  of  time  series 
that  are  plotted  in  Figures  12  through  17.  These  plots  were  generated 
with  pulse  inputs  of  length  T  =  {15,  20,  25}  and  W  =  -1.0.  The  reason  for 
two  plots  for  each  case  (T  =  15,20,25)  was  to  demonstrate  that  each  simulated 
intervention  will  result  in  slightly  different  signals,  with  basic  similarities. 

These  results  are  very  similar  to  Figure  3  and  Figure  18,  actual  MMW  data. 

A  closer  look  at  Figures  18  and  19  indicates  an  intervention  duration  of 
approximately  T  58  25.  A  comparison  of  the  actual  target  data  (Figure  18)  and 
the  simulated  data  (T  =  25)  (Figure  16)  shows  remarkable  agreement.  Figure  20 
is  another  simulated  intervention  of  this  dynamic  model  (9)  with  a  pulse  of 
length  T  =  25. 


V.  SUMMARY 

The  MMW  non- target  data  of  August  1978  was  analyzed  by  the  Box  and 
Jenkins  time  series  approach.  This  analysis  resulted  in  an  ARIMA  (1,0,1) 
model  with  consistent  estimated  parameters  for  both  <J>  and  0  as  well  as 
2 

a£  =  N(0,^  )  (white  noise).  The  ARIMA  (1,0,1)  model  characterizes  the  data 
remarkably  well. 

The  MMW  target  data  was  modeled  using  a  dynamic  intervention  approach. 
The  dynamic  model  proposed  is  promising  in  that  targets  can  be  characterized 
as  interventions  with  background  noise. 
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Figure  13. 


Simulated  Target  P^.  (15) 
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Figure  14.  Simulated  Target  (20) 
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Figure  15.  Simulated  Target  P.  (20) 
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Figure  16. 


Simulated  Target  P^(25) 
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Figure  17. 


Simulated  Target  P^. (25 ) 
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Figure  18. 


Actual  MMW  Target  at  30  Meters. 
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.  A  closer  Inspection  of  Figure  18. 
Actual  MMW  Target  at  30  Meters. 
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Figure  20. 


A  Simulated  Intervention  of  T  =  25. 
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